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Duringthe last decade, a great deal of attention has been paid to examining different quantum optical models 
[E S HO- [H. Recently an iteration technique @, @] has been suggested to solve the Schrodinger equation which 
Oh improves both analytical and numerical solutions of the problems and has been developed for some quantum optical 
systems. The solution of quantum optical Hamiltonians are important in the existing literature. In general, the study 
of two level-systems in one and two-dimensional geometry coupled to bosonic modes has been the subject of intense 
attention because of its extensive applicability in various fields of physics [1, 0, 11, 12, 13, 14, [3, 0, EH- 
qh Due to the practical and technological importance of these modelSjit is not surprising that various aspects have been 
studied both analytically and numerically III 2(1 21, 2^, 23 , 0, HE; 0, [23] • Such systems have often been analyzed 



' using numerical methods, because the implementation of the analytical techniques does not yield simple analytical 
expressions. Remarkably, exact solutions have not been thus far presented except for special cases, even though it has 
been suggested that the problem may be solved exactly. The required Analytical treatments need tedious calculations. 

Quantum optical models are one of the most fascinating phenomena in modern physics and chemistry, providing a 
general approach to understanding the properties of molecules, crystals and their origins. Most of these Hamiltonians 
have yet to be solved exactly. Therefore, the natural question arising at this point is: can the asymptotic iteration 
method be applied and used to obtain the solutions to these systems? The answer of this question is the main topic 
of this paper. 

In recent years much attention has been focused on the asymptotic iteration method (AIM) 0, [2^, 2jJ 3(J 31 , 3^] • 



This method reproduces exact solutions to many exactly solvable differential equations and these equations can be 
related to the Schrodinger equation. It also gives accurate results for non-solvable Schrodinger equations, such as 
the sextic oscillator, cubic oscillator, deformed Coulomb potential, etc. which are important in applications to many 
problems in physics. Encouraged by its satisfactory performance through comparison with other methods, we feel 
tempted to extend AIM to solve matrix differential equations. Although AIM has been applied to solve the Schrodinger 
equation, its application to the solution of the matrix equations [7| still needs to be improved. In contrast to the 
solution of the Schrodinger equation including potentials of Coulomb, Morse, harmonic oscillator etc. by using AIM, 
study of the matrix Hamiltonians has not attracted much attention in the literature. Therefore, we concentrate on 
the solution of a general two-dimensional two mode bosonic Hamiltonian by using AIM in this paper. 

The paper is organized as follows. In section 2 we discuss transformation of the Hamiltonians whose original forms 
are given as boson and fermion operators. We show that the Hamiltonian can be expressed as two coupled first 
order differential equations in the Bargmann-Fock space. In section 3, we develop AIM to obtain eigenvalues and 
eigenfunctions of the different matrix Hamiltonians widely used in physics. Section 4 is devoted to solve a wide class 
of physical Hamiltonians in the framework of the AIM. Finally we conclude our results in section 5. 
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TWO DIMENSIONAL TWO MODE BOSONIC HAMILTONIAN 

The most general form of the Hamiltonian consists of the coupling of a single spin-1/2 to the boson field in two 
dimensional geometry. It can be written as 

H = Ho + wocto + (kicl + n 2 a + + n 3 b + K4b + ) a+ + (71a + 72a + + 736 + 74& + ) c_ (1) 

where Hq = Tiui (a + a) + ?kc>2 (b + b) and Wj, Ki and 7j are physical constants. The Pauli matrices <ro.± are given by 

a ° = {o !) ; °"+ = (o o) ; ^ = (i 0)' (2) 

and bosonic annihilation a, 6 and creation a + , 6 + operators satisfy the usual commutation relations 

[a, a+] = [b, b+] = 1; [a, b] = [a+,b+] = [a+,b] = [a, b+] = 0. (3) 

It is seen that the Hamiltonian ([lj includes various physical Hamiltonians depending on the choice of the parameters, 
and it is Hermitian if k\ — 72, K2 = 71, K3 = 74 and K4 = 73 for real parameters. In general H describes spin phonon 
(photon) relaxation in the presence of a magnetic field. It can also be used to study interaction of a two- level atom with 
an electromagnetic field. As explicit examples, when lo\ = 102 = to, K2 = ftft = 71 = 74 = and K\ — K4 — 72 = 73 = k, 
the Hamiltonian H reduces to the E ® e Jahn- Teller (JT) Hamiltonian [3, |33j and when ui\ = ui + oj Ci u>2 = to — lo c , 
«2 = K3 = 71 = 74 = and K\ = — ka = 72 = —73 = k then the Hamiltonian H becomes the Hamiltonian of quantum 
dots including spin-orbit coupling One can also obtain Jaynes-Cummings (JC) Hamiltonian [3] and modified 

JC Hamiltonian [341 ] as well as many other interesting physical Hamiltonians by appropriate choices of the parameters 
u>i, Ki and 7^ in JT]). There exists a relatively large number of different approaches for the solution of the eigenvalue 
problem for Hip = EtJj in the literature. However, we present here a systematic treatment for the determination of 
the eigenvalues and eigenfunctions of |T]) in the context of AIM. 

We note that the Hamiltonian of a physical system is given in the form of a differential equation in some cases. One 
way to obtain the bosonic form of a Hamiltonian is to construct a suitable differential form of the bosons. Therefore, 
it is worth discussing some useful differential forms of the bosons and the connection between them before we begin 
to present a procedure to solve ([T]). 



Forms of bosons 



We are interested in the two-level systems in a one and two-dimensional geometry whose Hamiltonians are given in 
terms of bosons-fermions or matrix-differential equations. By the use of the differential form of the operators one can 
easily find the interrelation between boson- fermion and matrix differential equation formalisms of the Hamiltonians. 
There are various differential forms of the boson operators. At this point, let us start by introducing the following 
differential forms of the boson operators : 



2' 

b+ = i(x^ ///!--«'', -/■<>„)■ 1 J) 

'-!< 

where t = \/^P is the length parameter. In principle, if a Hamiltonian is expressed by boson operators, one could rely 
directly on the known formulae of the action of boson operators on a state with a defined number of particles without 
solving differential equations. Apart from the mentioned method, the Hamiltonians can often not be solved exactly 
then we need to develop alternative methods. It is amazingly interesting that the Hamiltonian ([T]) can be solved 
within the framework of the AIM when it is transformed into the form of the first order coupled differential equations 
under the constraints that will be given in Eqs. ()llaj) - (jlldj) . Now, we briefly discuss three well known differential 
forms of the bosons and the relation between them. In order to obtain different forms of the bosons we present a 
transformation procedure. For the sake of simplicity we take fi = m = lo = 1 and then £ becomes unity. 
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Consider the following exponential operator 

A = exp[/3(a+6 + &+«)] . (5) 

The operator acts on the bosons as follows: 

AaA' 1 = acos/3 - 6sin/3; Aa+A -1 = a + cos /3 - b + sin/3 (6) 
A6A -1 = b cos (3 + a sin (3; A6+A -1 = b + cos/3 + a+ sin/3 

If we set /3 = — 7r/4 and change the variable y — ► ij/ in the transformation given in ©, we obtain the following relations 

Aa+A- 1 = -L^-flS,); AaA" 1 = -^(s + d x ), 

Ab+A- 1 = j=(v-9y); AbA~ l = -j=(y - d y ). (7) 

The other important differential forms of the bosons can be obtained by transforming the bosons with the following 
operator: 

r = exp[|(a 2 + a+ 2 + 6 2 + 6+ 2 )] . (8) 

where a is the rotation angle of bosons. The action of the operator on the bosons is given by 

Tar -1 = a cos a — a + sin a; ra + T _1 = a + cos a + a sin a. 

r&r -1 = b cos a - b + sin a; rb+r -1 = 6+ cos a + &sina. (9) 
For a — 7r/4, the boson operators take the form 

a — > — — (a — a + ) = — — ; a + — > — — (a + a + ) = x 

The last formulation is known as the Bargmann-Fock space description of bosons [35| and this form plays a key role 
when constructing the one- variable first order matrix differential equation form of |T]). When we insert the form (|10[) 
in (fT]) , the Hamiltonian H becomes a first order and two variable matrix differential equation. In order to separate the 
variables, we look for the conserved quantity of the system. After some treatment, we obtain the following conserved 
quantities: 

(11a) 
(lib) 
(11c) 
(lid) 

We note here that although we obtain four conserved quantities, K\ conjugates with K 2 and N\ conjugates with N%, 
They can be provided by a similarity transformation 

K = K 2 = UK 1 U~ 1 and N = N 2 = UN 1 U~ 1 (12) 

where U = 0+ + o~_ . It is well known that if two quantum mechanical operators commute then they have common 
eigenfunctions. Thus one can write 

K\m n 2 ) = (k+ i J \m n 2 ) , N |m,n 2 ) = (k + - J |m n 2 ) , (13) 
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and the eigenvalue problem can easily be solved in the Bargmann-Fock space. Thus, we obtain the following expressions 
for the eigenfunctions of K and N: 



tp(x,y) 
ip(x,y) 



= x k 4>(xy) |t) +x k+1 <j){xy) |1) , for K 
= x k+1 cj)(y/x)\V+x k <b(y/x)\i) for N 



(14a) 
(14b) 



where ||) stands for the up state and ||) stands for the down state. In this case one can normalize the wavefunction 
tp(x,y) by using the relation 



-w 



ip(x, y)dxdy 

where W is the weight function and the integral is taken over all space. The eigenfunction of the Hamiltonian can be 
obtained from the relation 



|ni,n 2 ) =T 1 A l tp(x,y). 



(15) 



Meanwhile we note that the conserved quantity N crucially depends on the conservation of the number of particles. 
The classical motion of the particle takes place in the space of angular momentum on a sphere. However, to obtain 
a physical meaning of K from another perspective; when the motion of the particle takes place on an ellipsoid, 
the conserved quantity is given in the form of (|13p . We have obtained four different conserved quantities for the 
Hamiltonian (flj depending on the choice of parameters. However, it is fact that K\ and Ki can easily be mapped 
into each other. It implies that Hamiltonians obtained under the constraints (lllaj) and (lllc|) correspond to physically 
similar systems. We also say that iVi and N2 can also be mapped into each other and the Hamiltonians obtained 
under constraints (|llb[) and (jlld[) also correspond to physically similar systems. Therefore, it is valid to discuss the 
solution of the Hamiltonian (fTJ) under the conditions given in (| 1 lc|) and (jlld[) . 

Since Eqs. (jllap and ([1]) commute, then they have the same eigenfunctions under the constraint given in (|llap . 
Thus, insertion of (|14a|) into the Hamiltonian |T]) and using the form (fT0|) . we obtain the following set of one variable 
coupled differential equations: 



(UJI + L0 2 ) Z-J- + (k + i ) a.' I 



- ujq — E 



0i (*) 



d 

l + z Tz ]l - 



K4.Z 



02 (z) = 



(16a) 



(o>i +uj 2 )z — 
dz 



k + - ) wi 



id'2 — LjJq — E 



4>2 (z) + 



^ 3 dz 



■72 



01 (z) = 



(16b) 



where z = xy and E is the eigenvalue of the Hamiltonian H and <f>\ (z) and 02 (z) correspond to up and down 
eigenfunctions of the Hamiltonian H, respectively. Similarly, when we substitute (|14bj) into the Hamiltonian {T]) with 
the form (jTDJ), we obtain the following set of one variable coupled differential equations: 



,W2 



0J1) z— + (k + l)^i + ojq 
dz 



E 



01 (z) + [K2 + 4>2 (z) = 



(17a) 



(u>2 — Wl) Z— h ku>\ — LUq — E 

dz 



02 (z) + 



7i(fc + 1) + (73 - 7iz) — 
dz 



0i (z) = 0. 



(17b) 



where z = y/x. Our task is now to apply the AIM to solve the corresponding Hamiltonian. 



DEVELOPMENT OF THE AIM FOR MATRIX HAMILTONIANS 



In this section we systematically present a procedure for the solution of 2 x 2 first-order matrix differential equations. 
Consider the following first order matrix differential equation: 

0' = u cf> (18) 

where <fi = [0i, 02] T , two component column vector uq is a 2 x 2 matrix function. Note that and uq are functions 
of z and 0' is the first derivative with respect to z. Now, in order to obtain a general solution to this equation in the 
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framework of the AIM we use similar arguments to those given in Q. The differential equation (fT8|) can be written 
as the two coupled equation 

4>'i = «O01 + &O025 4>2 = C O02 + d o 0i (19) 

where ao, b , Cq and do are elements of the matrix uq. It is easy to show that the n th derivative of 0i and 02 can be 
written as 

0'/ = ai4>\ + bi4>2; 02 = c i^2 + di4>! 
4>" = a 2 4>\ + b 2 4>2; $2 = c 2 4>2 + d 2 4>i 

(20) 

cj)^ = a n -i4>i + b n -i4>2; 4> 2 n ^ — Cn-1^2 + dn-i0i 

= CL n <t>l + b n <f>2] 4>2 1+1) = C n4>2 + d„4>i. 

In order to discuss the asymptotic properties of ([1]), it is necessary to determine the coefficients a n ,b n ,c n and d n . 
After some straightforward calculation, we can obtain the following relations: 





— flO a n-l " 


r- a'n-1 ■ 


-1- do&n-l 
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= b()a n -i - 




- c 6 n _i 




= c c„_i H 




- 6o<in-l 


d„ 


= doCn-l - 


K4-1- 





(21) 

Our task is now to introduce the asymptotic aspect of the method. For this purpose, the n th and (n+ l)*' 1 derivatives 
of 0i and <f>2 can be written as 

A n ) (a , j. \ An) ( , , d n-i , 

01 = On-l 01 H 02 , 02 = Cn-l 02 H 01 

V a «-l / " V °n-l 

0i" +1) = a„ L + % 2 ) , 2 " +1) = c„ f 2 + ^0^ . (22) 

The coefficients do and Co include the coupling constants. Therefore, for sufficiently large n we can suggest the 
following asymptotic constraints: 

b n -\ b n d n -i d n 

— — — Ai; — — — A 2 . (Z3) 

In this formalism, the relations given in (|23[) imply that the wave functions 0i and 02 are truncated for sufficiently 
large n and the roots of the relations given in (|23[) belong to the spectrum of the matrix Hamiltonian. Therefore, 
one can easily compute the energy of the Hamiltonian by solving (|23[) for the energy term when z — > Zq. Under the 

,177771 • <h (n + 1) <i ( " +1) 

asymptotic condition of (|23p. one can find the wave functions 0i and 02 . When we take 1 M and 2 , , by using 

<t>\ <t> 2 

(122p under the constraints given in ((23|) . we obtain: 



exp I y ^~ dz J or ^2 = cx p ( J ~f— dz ) ■ ( 24 ) 

Substituting the expression of a n (and c„) given in (|2"Tj) into (f2~i| and then replacing the 1 ™' ) (and 2 ™' ) ) in (|20[) . 
respectively, one gets the following expressions; 



0i + Ai0 2 = exp I J (a + Xido) dzj or 2 + A 2 0i = exp I J (c + A 2 &o) dzj . (25) 
Using the second equality in (|25p . one can substitute 02 into the first equality in Eq. (|19[) . Thus, one writes 



0i + (b X 2 - a o )0i = Ci6 exp (J (c + A 2 6 ) dz 



(26) 
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and the solution is found as 



0! = exp (J (a - X 2 b ) dzj J C\b a e (f(co+^b )dz) dz + c ^ 
If the same procedure is performed for fa, one finds the solution as 

4> 2 = exp (J (c - Aidb) dzj J C 3 d e (f(^+Xtd )dz) dz + ^ 



(27) 



(28) 



An immediate practical consequence of these results is that the eigenvalues and eigenfunctions of the various 
quantum optical Hamiltonians can be determined. In the following sections, it is shown that this asymptotic approach 
opens the way to the treatment of a large class of matrix Hamiltonians of practical interest. 

RESULTS AND DISCUSSIONS 

In this part, we apply the results of previous sections to obtain the solutions of the Hamiltonians given by Eqs. 
(116al b) and (|17ai bV We briefly discuss the corresponding physical system of each Hamiltonian. 



Solution of the Hamiltonian H under the constraints: 71 = 74 = K2 = 13 = 

In this case the Hamiltonian includes two important physical Hamiltonians: E®e Jahn- Teller (JT) Hamiltonian § 
and Hamiltonians of quantum dots including spin-orbit coupling 1^, 33[. When u)% — u}% = u, K2 = K3 = 71 = 74 = 
and Ki = K4 = 72 = 73 = n, the Hamiltonian H reduces to the E <£> e JT Hamiltonian as we have mentioned before. 
It is obvious that the corresponding first order differential equations (|16a[ b) for JT Hamiltonian can be written in 
the following form: 



<t>'i = a a (f>i + b 4> 2 ; (f>' 2 = co02 + d (f>i. 



(29) 



where the coefficients ao, bo, cq and do are given by 



a Q = 

bo = 

co = 

do = 



k 2 -2uj -2k + 2E-2 

Az — K 2 
k (u! + k + 2z + E) 
k 2 - 4z 

k 2 (1 + k + z) + 2z (cj - k + E - 2) 

z (4z — k 2 ) 
n(E-Lj a - k + 2z + 1) 
z (k 2 — Az) 



(30) 



Using a simple MATHEMATICA program one can compute a n ,b n ,c n and d n by using the relations given in 
(f2T1) . On the other hand, for each iteration the expression 8\{z) = b n -i(z)a„(z) — a n -i(z)b n (z) (and 62(2) — 
d n -i(z)cn(z) — c n -i(z)d n (z)) depends on different variables, such as E n ,K,ujo and z. It is also noticed that the 
iterations should be terminated by imposing the quantization condition Si(z) = 0, i = 1,2 as an approximation to 
(|23f to obtain the eigenenergies. The calculated eigenenergies E n by means of this condition should, however, be 
independent of the choice of z. The choice of z is observed to be critical only to the speed of convergence of the 
eigenenergies, as well as for the stability of the process. In our study it has been observed that the optimal choice for 
z is z = |31J. Therefore, we set z = at the end of the iterations. We also note that the first value of the solution 
set of 5i(z) = (or 82 {z) = 0) is not physically acceptable unless the system is exactly solvable. 

To fix the iteration number for convergence, the first twenty energy levels have been determined for n = 
8, 9, 10, 11, 12, 14, 16, 18 iterations for the Hamiltonian above. It has been obtained that E20 = 21.103745 for n = 10, 
E 20 = 21.007 1 71 for n = 11, E 2 o = 21.007064 for n = 12, E 2 o = 21.007064 for n = 13, and E 20 = 21.007064 for 
n = 14, 16, 18 iterations, respectively, for k = 1, Wo = 0, w = 1, K = It is obviously seen that there is no change in 
the 20 t/l energy value for n > 12. Since it is also the same for the other Hamiltonians given below, n = 14 iteration 
is assumed to be sufficient for the determination of the energy eigenvalues of the related Hamiltonians. 
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TABLE I: The ground-state energies obtained from Refs. [3r3 . |37| . |3S| . [39l | and from the quantization condition by using 
coefficients given in Eq. (|30p for different values of k, in the case that k = 0, too = 0, it! = 1. 

k Ref. [36] Ref. [37] Ref. [38] Ref. [39] Present Results 



0.25 


0.774 


0.7766 


0.7765 


0.7739 


0.7738 


0.5 


0.578 


0.5877 


0.5870 


0.5780 


0.5780 


0.75 


0.400 


0.4173 


0.4158 


0.3998 


0.3997 


1 


0.233 


0.2586 


0.2560 


0.2331 


0.2330 


2 


-0.369 


-0.3157 


-0.3232 


-0.3686 


-0.3689 


3 


-0.919 


-0.8466 


-0.8575 


-0.9177 


-0.9189 


5 


-1.961 


-1.8716 


-1.8831 


-1.9540 


-1.9610 


7 


-2.976 


-2.8833 


-2.8932 


-2.9586 


-2.9760 


to 


-4.485 


-4.3937 


-4.4019 


-4.4594 


-4.4850 


15 


-6.991 


-6.9042 


-6.9108 


-6.9610 


-6.9901 


20 


-9.493 


-9.4111 


-9.4168 


-9.4627 


-9.4809 


30 


-14.496 


-14.4202 


-14.4249 


-14.4651 


-14.488 



a 

w 2 



k=0, coo=0.5 



0.5 1 1.5 

Coupling Constant, K 



FIG. 1: Energy of the displaced coupled harmonic oscillator as a function of coupling constant. 



In Table U we have compared our results with previous studies. It is seen that the results obtained by the AIM 

agree with those in Refs. [361 [37 , Sinn. 



Under the given conditions the Hamiltonian (p} takes the form: 

H = H a + uj <t + n [(a + 6+) cr+ + (a+ + b) aJ\ . (31) 

As a related topic, we mention here that studies of the E ® e JT problem led Judd f§| to discover a class of exact 
isolated solutions of the model. To determine the relations between the parameters of the model, one can obtain an 
analytic form of two eigenvectors of the Hamiltonian corresponding to the specific energy. The complete description 
of these solutions has been given by Koq et al. [3^|. They observed that quasi-exact solutions can be obtained by 
using osp(2, 2) super algebra. 

The physical systems described by the Hamiltonian (pTTj) are summarized as follows. When uiq — 1/2 and k = 0, 
the corresponding equation related to the Hamiltonian of the displaced coupled harmonic oscillator whose eigenvalues 
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k=0.5, coo=0.5 



0.5 1 1.5 
Coupling Constant, K 



>> 6 
& 4 



k=l, coo=0.5 



0.5 1 1.5 2 
Coupling Constant, K 



FIG. 2: Energy of the octahedral JT systems as a function of coupling constant. 



are obtained as a function of coupling constant, k, is solved by using AIM and the result is given in Figure 1. 

The Hamiltonian corresponds to three octahedral JT systems when uj n = 1/2 and k takes integer or half- integer 
values. These octahedral systems are r 8 ® (e + T2) linear E ® e and linear r 8 ® T2, for which their eigenvalues are 
depicted in Figure 2. 

Finally, in the presence of an external field, luo 7^ 1/2, the Hamiltonian is compatible with the generalized E ® e JT 
system and dimers. For dimers k = and for generalized E ® e JT system, fc takes half-integer values. The results of 
the iteration for dimers are given in Figure 3 and for generalized E ® e JT system are given in Figure 4. 

It is shown in the first 4 Figures that the energy becomes an oscillating function of the coupling constant n when 
k > 0. 

The Hamiltonians of quantum dots including Rashba coupling can be obtained when we set aj\ = LU2 = K2 = 
«3 = 71 = 74 = and Ki = — K4 = 72 = —73 = k, and the coefficients of the coupled differential equations are given 

by 

k 2 - 2w - 2k + 2E - 2 

Az + K 2 
k (2z - loq — k - E) 

k 2 +4z 
K 2 + 2 (lu + k + E) k + l 



k (E — oj — k — 2z — 1) 
z (k 2 + Az) 

Note that the origin of the Rashba spin-orbit coupling in quantum dots is due to the lack of inversion symmetry 
which causes a local electric field perpendicular to the plane of the heterostructure. In the literature, the Hamiltonian 
has been formalized in the coordinate-momentum space leading to a matrix differential equation. 



a = 

bo = 

co = 

d = 
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k=0, 0)o=0.0 



W 2 




>h 4 

u 





0.5 1 1.5 2 
Coupling Constant, K 

k=0, 0) = 1.0 



0.5 1 1.5 2 
Coupling Constant, K 

FIG. 3: Energy of the dimer as a function of coupling constant. 



Solution of the Hamiltonian H under the constraints: 71 = 73 = K2 = K4 = 

Under the constraints given in this section, we study various well-known exactly solvable Hamiltonians which give 
us opportunity to test our approach. By the given constraint, the Hamiltonian {1} takes the form 



H = H a + cj a + («ia + n 3 b) a + + (720+ + 746+) er_ 
We interpret below the solution of the three physical systems using AIM. 



(33) 



Jaynes Cummings Hamiltonian (74 = K3 = 0J2 ~ 0, K\ — 72 = k) 

The Jaynes-Cummings (JC) Hamiltonian with rotating wave approximation is given by 

H = uia + a + WqCq + k ((T + a + cr_a + ) 
In this case the coefficients of the coupled differential equations (|17al b) are given by 

' k 2 + lu (E — kuj + u>o) " 

ao - 
b 



k{E- uj a ) 



Co 



LU Z Z 

—E + ll> + ku> + luq 



(34) 



(35) 
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FIG. 4: Energy of the generalized E ® e JT system as a function of coupling constant. 



The AIM leads to the following expressions for the eigenvalues of the JC Hamiltonian 



n = 1 

n = 2 
n = 3 



E = 

£ : 

£ : 
E 



1 



1 



* + - J w T ^V 4k2 ( fc + !) + ( w + 2cJ o)' 



k--jiuT g V 4 ^ 2 (* - !) + (w + 2^o) 
^fc + | - w =p ^ (fc + 2 - n) + + 2w ) 2 



(36) 



It is obvious that when the coupling constant k is zero, then the result is the eigenvalues of the simple harmonic 
oscillator. 

Substituting the variables oq, bo, Cq and do in Eqs. (f2"T)l and 



n 
n 
/? 



1; X = 1 

2; 0i = z 



one finds the eigcnfunctions 4>i(z) and 02 (^) a s 
2 = 1 



= 3; 



Using Eq. (fl"4b"|) , one writes 



and finds 



n = n; <f> 1 = z™" 1 cf> 2 = z n ~ x . 
iP(x 1 y)=x\y/xT- 1 \l)+x k + 1 (y/xT- 1 \ T) 
iP(x,y)=Cx 1+k - n y n - 1 (\ i)+x\ t» 



(37) 
(38) 
(39) 



where C is the normalization constant, k is the state number and n is the iteration number. In order to find the 
original eigenfunction of the Hamiltonian, one can use Eq. (|T5]) by replacing the operators given in (0) and ©■ 
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FIG. 5: Energy of the Rashba Hamiltonian as a function of coupling constant. 



Dirac Oscillator(uj\ — ui2 — ^3 = 74 = 0) 

One can also show that the constrains given in this section include the Dirac oscillator. In order to show this we 
express the Dirac oscillator with boson operators. Consider the (2 + 1) dimensional Dirac equation for a free particle 
of mass m in terms of two component spinors, then ip can be written as 



Etp = m iVi + mc 2 c7 j ip. (40) 

The momentum operator pi is a differential operator p = —ih(d x , d y ) and the 2D Dirac oscillator can be constructed 
by changing the momentum as p —> p — imcoaor. Then the Dirac equation (|40p takes the form 



(E — mc 2 ao) = c [(p x — ip y ) — imuj'(x — iy)] <j + + c \(p x + ip y ) — imuj'(x + iy)\ cr_ . (41) 
After some straightforward treatment, we obtain the bosonic form of the Dirac oscillator: 



(E — mc 2 (7o) ip — 2icV muj'h [aa+ + a + cr_] -0. (42) 

An immediate practical consequence of these results is the Lie algebraic structure of the Hamiltonians that can easily 
be determined. It can be obtained by setting the parameters of the JC Hamiltonian to uo = 0, k = 2icV rruv'h and 
luq = mc 2 . Then the eigenvalues of the Dirac oscillator are given by 

E = ±^4m 2 c i - Ahco'mc 2 (k ± nj. (43) 



Modified Jaynes-Cummings Hamiltonian(K,i = 72 = K3 = 74 = K, ui2 = u\ — lu) 



In addition to these Hamiltonians, we can also show that our formalism includes another important Hamiltonian: 
When a single two-level atom is placed in the common domain of two cavities interacting with two quantized modes, 
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the Hamiltonian of such a system can be obtained from the modification of the JC Hamiltonian and is given by 

H = u;a + a + ujb + b + lu ctq + k (a + b) a + + k (a + + b + ) cr_ . (44) 
Without detailed calculations, one can obtain the energy of the MJC Hamiltonian in the closed form by using AIM: 

E= + ± + l -n)n 2 + {2uj - l) 2 . (45) 
Consequently the exact eigenvalues can be reproduced by AIM. 

CONCLUSION 

We have systematically discussed the solutions of various physical Hamiltonians within the framework of AIM. We 
have shown that the formalism given in this paper leads to the exact or approximate solution of the problems of various 
physical systems. We have applied the AIM to the problem of an electron in a quantum dot in the presence of both 
a magnetic field and spin-orbit coupling. The procedure presented here gives an accurate result for the eigenvalues of 
the both JT and Rashba Hamiltonians. The suggested approach can easily be used to solve other quantum optical 
problems which are not discussed here. 

We have presented a transformation procedure that offers several advantages, especially if one wishes to describe 
the eigenvalues of the bosonic Hamiltonians by using AIM. We have also presented the steps towards an extension of 
the AIM. 

The technique given in this article can be extended in several ways. The Hamiltonian of a quantum dot including 
position dependent effective mass may be formulated and solved within the procedure given here. We hope that our 
method leads to interesting results on the spin-orbit effects in quantum dots in future studies. Along this line we have 
work in progress. 
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